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ABSTRACT 


The authors have developed an interactive graphics system for 
the correction and reduction of airborne laser profiles of sea ice. 
The system is implemented on a Vector General renewable display system 
connected to a PDP11/45 computer, but minor modifications would make 
it suitable for any computer with interactive graphics. The laser 
profile is displayed in sections and the human operator is armed with 
routines for the removal of noise spikes, phase jumps and aircraft 
altitude variations. Completely automatic profile correction techniques 
have been found to be inadequate or unworkable, especially in heavily 
ridged ice. Although designed specifically to treat the output of the 
Geodolite laser profilometer, the system can be applied to many other 
data correction problems in geophysics and other sciences where long 
time series of data are obtained with periodic noise spikes or a long- 
period drift or modulation which requires removal. 
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1. INTRODUCTION 


A knowledge of the surface topography of sea ice is of great import- 
ance in the modelling of polar ocean dynamics. The drag coefficient 
between wind and ice is determined by the surface roughness (Banke et 
al, 1976), and since most of the heat flow from ocean to atmosphere 
occurs through thin ice (Maykut, 1976) it is valuable to measure the 
distribution of surface elevations as a guide to the distribution of 
ice thicknesses. The airborne laser profilometer is the most effective 
instrument for swiftly generating such data over long stretches of the 
ice cover. It produces a profile of surface elevation along a track 
directly below the aircraft, using a beam whose footprint is essentially 
a moving point and which is therefore capable of very high resolution. 


The earliest and most extensive laser surveys of the Arctic Ocean 
were carried out by the US Naval Oceanographic Office (Ketchum, 1971; 
Welsh and Tucker, 1971; Tooma and Tucker, 1973) and furnished data 
which were used to test theories of pressure ridge height and spacing 
distributions (Hibler et al, 1972, 1974; Diachok, 1975) and to delineate 
regions of the Arctic Ocean with differing degrees of roughness (Tucker 
and Westhall, 1973; Hibler, 1975). The same agency has carried out 
laser profiling in the Beaufort Sea as a contribution to the AIDJEX 
project, and co-operative profiling with a submarine in the East Greenland 
current (Kozo and Diachok, 1973). Laser profiling in the Arctic Ocean 
has also been carried out by an Argus aircraft of Department of National 
Defence, Ottawa (Dunbar and Lowry, 1974; Wadhams, 1976). As well as 
measuring ridging, it is also possible to use a laser profilometer to 
measure attenuation rates of ocean waves in the marginal ice zone 
(Wadhams, 1975). To date no laser profiling has been done in the 
Antarctic, although sea ice characteristics there differ significantly 
from those in the Arctic and a profiling experiment would be of great 
interest. The most recent laser profiling experiment in the Arctic 
Ocean has been a co-operative project in October 1976 involving an 
Argus aircraft and a British nuclear submarine following identical 
tracks at approximately the same time (Wadhams, 1977; Wadhams and 
Lowry, 1977), thus enabling the statistical relationship between surface 
and bottom topography to be derived. The interactive graphics system 
described in the present report was developed to analyse the data from 
this experiment. 


So far all airborne laser surveys have employed the same instrument, 
the Spectra-Physics Geodolite 3A. This uses a continuous He-Ne laser 
beam with a wavelength of 632.8 nm and an output power of approximately 
25 mW. The outgoing beam is sinusoidally modulated by an electro- 
optical modulator, and after reflection from the target the beam is 
received by a Schmidt-Cassegrain telescope of 0.2 m diameter which is 
coaxial with the laser. The phase shift between the modulations on 
the transmitted and received beams is measured, giving the range to the 
target. A set of modulation frequencies is available, usually differing 
by factors of 10. This enables the laser to be used as an accurate 
altimeter: the lowest modulation frequency gives a full scale deflection 
(O to 27 in phase) for an altitude change of 10° feet, and the higher 
frequencies give 104, 103, 102 and 10 feet full scale deflection. 


By switching frequencies stepwise, each significant figure of the altitude 
in feet can be obtained. When used to obtain sea ice terrain profiles 
only the higher frequencies are used, giving better resolution. 
Unfortunately, the aircraft normally porpoises with an amplitude exceeding 
10 feet, and ice ridges can also exceed 10 feet in elevation, so that 

a laser profile on the 10 foot scale has frequent "phase jumps". A 

phase jump occurs when the range crosses one end of the 10 ft interval, 

ie the phase crosses 0° or 360°. When this happens the electronics 
automatically put ina 360° jump, positive or negative, in order to put 
the profile back on scale. 


Figure 1 (top) shows a typical raw profile produced by this type 
of instrument. There are four components in the profile:- 


1. The actual sea ice elevation. 

2. A profile of aircraft porpoising, in which both altitude and 
attitude variations combine to give a smoothly-varying range oscillat- 
ion of wavelength about 2-4 km. (ie Figure 1, centre). 

3. Phase jumps, which are theoretically instantaneous and of identical 
amplitude. In practice, because of the recording system used 
(magnetic tape) and imperfections in the profilometer, the jumps 

have a finite rise time of up to 0.05 s and are of varying amplitude. 
The profilometer used for the 1976 experiment was modified to give 

a 20 ft full scale deflection on its most sensitive scale, reducing 
the frequency of phase jumps, and to give a jump of 180° instead of 
360°. This restores the profile to the centre of the scale instead 
of to the opposite extremity, reducing the likelihood of the profile 
immediately jumping back by the same 360°. 

4. Noise spikes and "drop-outs" (none visible on fig. 1). A drop- 
out is a brief positive-and-negative spike caused by a momentary loss 
of phase lock. The loss may be due to the aircraft crossing an ice- 
water boundary, or else experiencing an exceptionally sudden range 
change (eg an overhanging ice block in a ridge). Other spikes may 

be caused by interference from other equipment on the aircraft. Longer 
periods of noise can occur due to phase lock failure over thin 

cloud or undercast mist. 


The purpose of a data processing system is first to remove (3) 
and (4), leaving a profile as in fig. 1 (centre), and then to remove 
aircraft porpoising as accurately as possible so as to leave a true 
profile of ice surface elevation as in fig. 1 (bottom). Computer 
programs can be written to automatically remove phase jumps, but the 
process is very unsatisfactory since the profiles must always be 
examined afterwards and sections with imperfect jump removal or with 
noise spikes present must be discarded or further treated. Similarly, 
there is no fully satisfactory way of automatically removing aircraft 
motion. A simple filter will not work, since the data do not represent 
a gaussian random process, the deviations from the porpoising being all 
in the positive direction. A three-stage filter was described by 
Hibler (1972) to overcome the last problem. The method used is to connect 
up the local minima, thus underlining the data. This underlining is 
then filtered and subtracted from the corrected data to yield an ice 
profile. The main disadvantages of this procedure are that it does not 
allow for the existence of leads, and again that it does not allow for interactive 
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A laser profile of sea ice, 
(c) with aircraft porpoising removed. 
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viewing. Thus it must be coupled with an automatic phase jump/noise 
spike removal program, which works only in regions of good quality data 
(ie, normally lightly ridged areas). 


The system described in this report is based on that developed at 
Defence Research Establishment Ottawa by Lowry and Brochu (1976) and 
uses a human operator to do the pattern recognition for both errors and 
aircraft porpoising removal, leaving the computer to manipulate the 
data and do all the necessary computations. Data is presented on a 
television-type display and the operator indicates the position and 
nature of errors with a light pen. When the errors are removed to the 
operator's satisfaction the data is underlined by the operator, again 
using a light pen. The underlined data is then filtered and subtracted 
from the corrected data, leaving a profile of pure ice elevation, free of 
aircraft motion. The corrected profile is written to magnetic tape so 
as to be available for statistical analysis on the same computer or on 
a separate general-purpose computer. 


2. DESCRIPTION OF THE PROGRAM 


The system has been implemented on a PDP11/45 computer, operating 
under DOS, with a disc drive, two magnetic tape drives and a Vector 
General visual display unit with light pen and keyboard. The program 
was written mainly in Fortran. 


The program sequentially displays blocks of data from an input 
tape; these are interactively corrected by the operator and the 
result displayed for inspection. If satisfactory, the blocks are 
output onto a second magnetic tape. During this process, sufficient 
data is kept on disc so that the operator can stop at any time and pick 
up later where he left off, and so that there is minimum loss should 
there be a system crash. 


The main program keeps track of the block numbers, controls the 
disc and tape files, and calls the correction subroutines. After 
asking the operator whether the data is new or old (ie system being 
restarted after a break) the program reads in the relevant block of 
data into the array ICE, which also holds subsequent corrected versions 
of the data. The program then hands this over to the subroutine 
CORREK for the correction of phase jumps and noise spikes, and then to 
the subroutine UNDERL for the drawing of the aircraft motion. This 
second routine produces an array PTS holding the drawn underline. 


This underline, which consists of line segments, is filtered by 

a Hamming filter, which essentially produces a weighted running mean 
of 200 points. However, to avoid edge effects of the filter at the 
beginning and end of a block, the array PTS must be preceded by at 
least 100 points of the underline from the previous block (stored in 
array BSTR), and must be followed by the first 100 points of the subsequent 
block (stored in array ESTR). The very first and last blocks on the 
input tape must therefore have end ramps constructed for them. For 
this reason the filtering (subroutine HAMFIL) must take place a block 
behind the correction and underlining. Thus for new data the sequence 
of events is:- 

-- process blocks 1 and 2 

-- filter and display block 1 

-- process block 3 

-- filter and display block 2 etc. 


When a block is filtered and displayed, the operator may reject it, 
in which case this block is treated as a first block of data and the 
program restarted, or he may accept it, in which case the smoothed 
underline is subtracted and the result similarly displayed for approval. 
When finally approved, this block is written onto output tape. 


It should be noted that the program has been linked with three 
standard libraries of routines held permanently on the Cambridge 
PDP11/45 computer:- 


1. A standard Fortran subroutine library 
2. Magnetic tape routines 


MTREAD reads a block 
MTWRT writes a block 
MTWIM writes a tape mark 
MTSKIP skips blocks 
MTRWND rewinds 
3. Graphics routines 
NITDEV initialise device 
NITBUF initialise a picture buffer 
Saas | for creating a picture segment 
DELPIC for deleting a picture segment 
WRITOL for positioning tracking square 
INWAIT for inputting on visual display unit/Vector General 
keyboard 
REATOL for finding position of tracking square 
5 ae absolute and relative line 
LINIR segment drawing 
WINDW for setting up a display window 


The subroutine CORREK has two main functions. Firstly it seeks 
large jumps between consecutive points. It indicates these in turn, 
asking the operator to say ':- 


i) that it is a phase jump or noise spike, in which case it is removed, or 
ii) that the jump is genuine data, in which case it is kept. 
After this the routine waits for commands from the operator to remove 
remaining errors, redraw, reject and start again, or finally accept the 
modifications (see OPERATOR'S GUIDE). 


The subroutine UNDERL waits for commands to insert or delete a 
line segment, alter the starting point of the line, or finally accept 
the underline. 


A basic block diagram for the main program is shown in Fig 2. 
Block diagrams for the subroutines will be found in section 4. 


Some modifications to the simplest form of the program were necessary 
because of space and speed limitations on the PDP11/45. The program fills 
most of the available 45K of store, so space saving procedures were 
employed which makes the program (section 4) slightly difficult to 
follow. Further HAMFIL as written in FORTRAN takes 90 secs for a filter, 
and so was replaced by an equivalent routine written in Assembler 
code, which requires only 15 secs. 


The Hamming Filter was chosen for its low side lobes in the 
Fourier transform domain. It should be pointed out however, that a 
simple top hat filter could be used. Its first side lobe is only 
13dB down, but in many cases this is entirely adequate. (If the data 
are to be subject to Fourier analysis, however, the side lobes become 
more important). A top hat filter, once filled, will filter any 
length of data block, with only a few simple calculations, thus taking 
much less time. 


3. AN OPERATOR'S GUIDE TO THE SYSTEM 


This section refers specifically to the operating system of the 
PDP11/45 - Vector General. It is included to illustrate the relative 
ease of carrying out these operations on an interactive graphics system. 
The program is set up on the user's disc; the data is in digital form 
on magnetic tape; the output is written to a second magnetic tape. 
Normally Geodolite data is recorded on analogue magnetic tape, so that 
A-to-D conversion is necessary. The 1976 data was recorded thus on a 
Tandberg Series 100 FM instrumentation tape recorder with flutter compen- 
sation, and digitally converted using an ADO1-D analogue-to-digital 
converter built into the PDP11/45 computer. Digital data points were 
0.01 s apart, corresponding to about 0.9 m of track; an analogue 
low-pass filter was included in the A-to-D converter to remove all 
energy above the Nyquist frequency. 
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Figure 2. Simplified block diagram of the main program. 


Initialisation 


The operator loads the input tape (data tape) on unit O and the 
output tape on unit 1. He loads the disc on the disc drive and initialises 
the PDP11l with the sequence 


HALT 
LOAD ADDRESS 173100 
LOAD ADDRESS 177406 
START 


The Vector General display is switched on and adjusted for high gain 
(since the light pen responds more readily at high intensity). 


When the PDP1l is initialised the operator logs in on the Teletype 
keyboard using his password (eg LO 61,166) and assigns the Vector 
General to output stream 8 


AS vG:,8J 
then starts the program running with 
RU G) 


where G is the name of the program. 


Operations 


1. Upon initialisation of the program the operator is presented with 
the question 


NEW DATA OR OLD? (N OR O) 


in the top left hand corner of the screen (where all comments or questions 
from the program to the operator appear). He answers on the Vector 
General keyboard with 


O) if he wishes to continue with a dataset begun in an 
earlier session, or 
NJ if he wishes to start a new dataset. 


For new data the program sets up its own ramp function to fill up one 
end of the Hamming filter. Otherwise it searches for the last block to 
be treated and uses the last part of that block to fill up the filter. 


2. The program now presents a block of 2000 points (about 2 km of data) 
on the screen, scaled so as to fill the active area, with a y=O line 
included (although this may be in an arbitrary position due to wrap- 
around) and with a vertical bar in the bottom left corner representing 
the scale of 200 digitising units (about 0.86 m). If the profile has 
any line segments for which 


y(n) - y(n-1) 3 200 


the tracking square will position itself at the bottom of the screen 
directly below the jump in question, starting with the jump nearest the 
start of the block. The operator now has the option of accepting the 
program's judgement, ie of having the jump removed on the grounds that it 
is a phase jump or a noise spike. To do this he just types 


) 


If he thinks, however, that the jump represents real data (part of a 
ridge) or that it can be better removed by. the "G" facility (see 3 below) 
he can reject the program's decision by typing 


N 


The tracking square now moves on to the next jump until all jumps are 
exhausted. An instantaneous phase jump has an amplitude of 10 feet 
(712 digitising units); we use 200 as the criterion because jumps 
seldom occur instantaneously. 


3. If there are now critical jumps in the block, or when the operator 
has dealt with all those that exist, the program presents on the screen 


TRY. GoR xX OR! © nnn 


where nnn is the number of the block (for information purposes}. If the 
Operator types 


RQ 


he can see the results of his operations in stage (2) since this redraws 
the block in its updated state, ie with the results of all manipulations 
to date. R can be used at any point during this stage of processing to 
examine the latest state of the data. 


The "G" facility deals with any undesirable features of the data that 
still remain. It replaces a short length of data with a horizontal 
straight line segment. It can therefore be used to draw a line to 
replace a noise spike, or else to take out a phase jump which occupies 
a finite number of data points instead of being instantaneous. To use 
this facility the operator uses the light pen to move the tracking 
square to the left side of the proposed line segment and types 


G) 
He then moves to the right side of the proposed segment and types 


2 


For typical applications, see EXAMPLES. 


If the operator is not satisfied with the alterations that he has 
made he can type 


X) 


which redraws the block in its original unmodified state, rejecting all 
of the operator's manipulations to date. 


When the operator is finally satisfied that the block now consists 
only of real data, with spikes and jumps removed, and subject only to 
aircraft motion (which should appear as a smoothly trending curve) he 


types 
oO. 


The program now proceeds to the next stage. 


4, It redraws the block once more and presents the message 
UNDERLINE 


It also presents the tracking square at the extreme left side of the 
screen on the y-value that it had reached while underlining the preceding 
block. If the operator is dissatisfied with this as a starting position 
for his underlining, he moves the square to a better position at the 
start of the block and types 


B) 
Otherwise he simply moves the square along (using the light pen), stopping 
at points which he feels are representative of the “level ice surface" 
and typing 


In this way he constructs a piecewise straight line profile of the surface 
of level ice. To recognise this surface requires judgement, experience 
and sometimes imagination; a prior study of some vertical aerial photographs 
of the Arctic ice cover is an advantage. The main problem is not to be 
misled by very smooth ice (with a "roughness" which is only the noise 
level of the original tape recorder). This is often thin polynya ice, 

and the operator should try to keep a short distance above such ice, 
following the level of smooth floe ice (which usually has greater 
roughness). Narrow cracks can also be mistaken for level ice segments. 
The aim is to achieve consistency; ie to follow either the level of 
polynya ice throughout (for open icefields) or to follow the level of 
smooth thick ice throughout, which is more appropriate for the central 
Arctic where polynyas may be far apart. 


If the operator is dissatisfied with his latest line segment he 
can type 


Dy) 


which deletes the last segment drawn. Repeated use will delete part or 
whole of his line. If by mistake he does not proceed monotonically in 
the positive x-direction he will receive the message 


MOVE ALONG 


and his last segment will be disregarded. If he exceeds the maximum 
allowable number of segments in a line (25) the program sends 


DELETE SOME 
and the operator must delete some of his points and redraw the rest of 


the line with points further apart. 


After taking his line to the extreme right of the block (with the 
help of a chart recorder printout, if available, to determine the 
continuing trend), the operator types 


oO) 


which sets the filter routine in action. 
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5. The program now filters the preceding block's underline, with the 
aid of the first few points in the current block to fill up the Hamming 
filter. When the filtering is complete (about 15 seconds) the program 
presents the preceding block together with its piecewise line filtered 
into a smooth curve, and asks 


OK? Y OR N 


The operator examines the smoothed underlining to see whether it appears 
to be a true representation of the level ice surface. If he thinks that 
it is, he types 


Y.) 
otherwise NQ 


A rejection takes the program back to stage (4) for the preceding block; 
it presents this block again for underlining. The formerly current 
block is lost and must be processed again after the preceding block. 

On acceptance, the program subtracts the smoothed piecewise profile 
from the laser data to obtain a final laser profile, allegedly a true 
representation of the. ice surface. It presents this again to the 
operator with the message 


OK? (Y OR N) 


Again the operator has the option of rejecting it, with the same consequences 
as before. If he accepts the final version, it is written onto the output 
tape and the program asks 


A REST? (Y OR N) 
If the operator wishes to wind up the session he types 
YQ 


and the program rewinds both tapes, having stored the current operating 
data. The operator now types 


FI) ; 
on the Teletype keyboard, which stops the program so that it is safe for 
him to remove tapes and disc. 
If he wishes to process another block, he types 
N) on the Vector General keyboard 


whereupon the next block of raw data is read in from the input tape and 
presented to him. 


il 


4, ILLUSTRATED EXAMPLES OF THE USE OF THE SYSTEM 


In this section we show several photographs of the Vector General 
display screen which demonstrate aspects of the operation. Figure 3 
shows the calibration signal used to eliminate errors from the tape 
recorder and digitisation. The profilometer output could be held at 
any one of three DC values corresponding to +180 offset, zero, and 
-180°. At least 5 seconds of each of the three signals were recorded 
at the beginning of each tape and whenever adjustments were made to the 
sensors. The left hand side of the trace shown in fig 3 has a 0° to -180° 
calibration signal. The average value of each offset was calculated 
from a print-out of the screen display, and used to establish the 
relationship between the unit of digitisation in the computer and the 
actual elevation perceived by the profilometer (the absolute accuracy 
of the profilometer is dependent on the stability of a quartz oscillator 
kept in a temperature controlled oven, and is of the order of 1 part 
in 10°). Figure 3 also gives an impression of the noise level introduced 
into the data by the recorder and digitiser; this is just the thickness of 
the "Line" representing the DC signal. 


Figure 4 shows a typical section of heavily ridged ice from the 1976 
co-operative profiling experiment. The vertical bar at the lower left 
corner is included by the program to give the operator an idea of the 
vertical scale, since the computer always re-formats the profile to 
occupy the whole height of the screen. The bar is 200 units (0.86 m) 
high. The tracking square has placed itself automatically beneath a 
spot where there is a jump of more than 200 units. In this case it has 
clearly identified the position of a 180° phase shift. 


In figure 5 the operator has accepted the program's judgement and 
by typing ) has removed the jump. The jump was a simple one, and its 
removal reveals the presence of a high pressure ridge. 


There is a second phase jump on the far right of the section. This 
has a finite rise time and was not found by the program. The operator 
moves the tracking square to the beginning and end of the jump, and by 
using the "G" command reduces the jump to a short horizontal line segment. 
Figure 6 shows the result of this operation. 


Since the raw data have no further phase jumps or spikes the operator 
now underlines it with short line segments joining apparent minima or 
stretches of smooth ice (figure 7). Figure 8 shows the result of filtering 
this piecewise line into a smooth curve. The operator accepts the result 
and subtracts the smoothed line from the profile. Figure 9 shows the 
final result. 
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Figure 3 


Figure 4 
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Figure 5 
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5. LISTING OF THE PROGRAM 


The program listing is followed by block diagrams of the two most 
important subroutines (figures 10 and 11); a block diagram of the main 
program was given as fig 2. 


15 


CORR ROK OO OOK OK IK OK IR OR OK I IO I I I IK aK I kK aca IC 2k 2 2k 2k 2k 2k 2k ak 2 a 2k ac aie ake ok ic ok 2 
(Ck KORO KOK KOK RIO OK I OK OK I OR IO I OK OK KK OK IK 2 OK kK 2k aK ak aK ok ak ak 2 ak 2c 2k 2k ok 2 aie ie a aie kc 2k ie 2c 2c ok 2 2c ac 
C MAIN PROGRAM FOR THE CORRECTION OF LASER PROFILES OF SEA ICE 
KOK KOK KOK KOK KOK OK OK OR OK IK OK I OK KK IKK OK aK kK kk I aK 2 aR aK aK a IC 2 2k a ak 2 2 2 aK 2k 2k ic 2k 2k 
COMMON /A/WC( 4) /B/TCC2) 
INTEGER ICEC 2000). PTS( 2000)» PICBUF(C 4200)s, BSTRC 101), ESTRC 101) 
INTEGER ANCIL(105)5Vs OFSET 
EQUIVALENCE (Ns ANCIL(1))5COFSETs ANCIL( 2) )s (LPTs ANCIL( 3) ) 
EGUIVALENCE(Ms ANCIL( 4) ).(BSTRC1)sANCIL(5))4¢(BSTRC 101), PTSC1)) 
EQUI VALENCECPTS( 2000). ESTRC1)) 
CK OK ROOK ORO IO OOK OI OK IKK OK IK IK IK IK I I IK a a I ko a kK aK I ak ak kK a aK ak aK ak aK ak ak ak ak 2k ak 2k oc 2k 2k 
E INI TLABIZALLON 
CK 2K 2K 2K 2 KK 2K 2K 2 2K 2K 2K 2K 2K 2 2K 2K oie 2k 2K 2k 2K 2 2K 2K fc ake 2 i 2c 2c 2 ik 2K 2K 2k 2K 2K IK 2k 2K 2K KK 2K 2K 2 2K 9K 2K 2K ic 2k 2K 2K KK 2c 2 2K 2c 2k 
DEFINE FILE 1¢€145105sUsNAN) 
DEFINE FILE 2¢€ 2, 2000s Us NPTS) 
DEFINE FILE 3¢ 2s 2000, Us NICE) 
CALL NITDEV(S) 
CALL NITBUFCPICBUFs 4200) 
CALL MTRWND(O) 
CALL MTRWNDO1) 
CALL HAMTAB 
CK KK 3K 2K CO IK 2K 2c 2 2 2K 2k ic 2 2K 2K 2K 2g 2 2K 2K 2K 2 2K 2 2K 2 2K 2 26 2K OK KK 2K 2K 21k 2c 2K 2 2K KK 2k 2k 2 2K 2k 2K 2k 2k 2K 2 2 2K 3K 2K OK 2K 2K 2K 2k ok ok ok ok 
C STARTING AFRESH OR CONTINUING ? 
C 2K 2K 2K 2 2g 2 2 OK 2K I 2 2 OK 2K 2K 2K 2K OK I 2K 2k 2K i fc ic ic 2c 2k 2K 2K 2 2G 2 2 2K IK 2g 2k 2k 2k 2k KK 2K 2K 2k ofc ot KK ik 2c 2K 2c 2K 2K 2K 2K 2k 2K 2 OK ok OK 2K 2K 
M=0 
1 CALL BGNPIC(C3) 
WRITEC&s 100) 
100 FORMATCIH1s,'"NEW DATA OR OLD? (N OR O)') 
CALE SENEPHEe 
READ(8, 101) V 
101 FORMATCA1) 
CALL ADELE RG@s) 
PEGs EQet Nip meGOLOme 
C2 2K OK 2K KK 2 OK 2c KK 2K 2K 2 2K OK 2 2 Ic 2K KK 2k 2 2K 2K 2K 2K 2k 2K oi 2c 2k 2K I 2g 2K 2k 2g 2K 2K 2g 2 og 2k 2k 2c 2K 2K 2k 2k a 2k 2k 2 aK 2k ok ok 2k 2k 2k ok ok ok ak ok 
C CONTINUING 
CK 2K 2K 2K OK 2 2K 2K 2K 2 2 2k 2 2k 2k 2 2 2k aK 2 2g 2 2 2c ak 2k i kc i 2 I 2 2k 2 2k 2 ako 2k ok 2k ak 2c ofc ok 266K 2c ok ak 2k ak ok 2K a 2k 2 aK 2k 2k 2k 2k ok ok ok 
READG!* 1S @ANCIE 
CALL MTSKIPC1sMs IERRs 1EOF) 
GOTO 3 
C2 2K 2K 2 2K 2K 2K 2 oc ic 2h 2 ic 3 ke 2K 2k 2 2 2c 2 IK 2c 2K Ic kg 2 2k ofc 2c kK 2c 2K 2K kc 2k 2 iC 2 oC 2k 2c 2K 2c 2 2k 2c ak akc 2c 2 2g 2K 2K 2K 2 a oi i ok ok ok ok 2k 2k 2k ok ok 
Ec STARTING AFRESH 
CK 2K 2K 2K OK 2 2K KK 2K 2K 2 2 go oc 2K oI 2 ok 2k 2K kc ic 2 oi ico 2c Ko 2K KC I IKK ok akc oc i aK ok 2k 2K 2c oe 2K ofc 2K 2K 2K 2K ofc 2k ok ok oo oe aK ok ok aK 2k ok 2k 
2 CALL BGNPICC3) 
WRITEC&s 106) 
106 FORMATC1H1s'STARTING AT BLOCK ') 
CALL ENDPIC 
READ(8s5105) N 
105 FORMATC( 14) 
CALL MDELPIeCGs) 
OF SET=0 
Beran 
oa CALL TAPINCICE,N»s1EOF) 
CALL CORREKCICEs OFSETsN) 
CALL UNDERL(PTSs LPT) 
CALL HAMBGN(BSTR) 
K=MOD(Ns 2) +1 
WRITECL 1 eANGLE 
WRITEC2'K) PTS 
WRITEC3*K) ICE 
CK ROOK ROOK KOO KOO I OKI OK KOK I I IK I KK IK I I I aK aK ak aK a aca ak ak ak 
C PROCESS NEXT BLOCK 
CK > 2K 2 2K KK 2 2K ko a I IK OK KI I I IK I A IK aI KOK KK II I RK KK RK KR oak fe a OK ok oko 
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N=N+1 

CALL TAPINCICEsNs I EOF) 
IFCIEOFe-EQ-W) GOTO 4 
CALL HAMENDCESTR) 

GOTO 5 

DO 41 K=1., 2000 
ICECK)=ICECK)+OFSET 

CALL CORREKCICEs OFSET,N) 
CALL UNDERLC(PTSsLPT) 
K=MODCN, 2) +1 

WRETEGS]* Ke PTs 
WHERE GS" hoa lL. CE 

CALL HAMCOPCPTS,s ESTR) 
K=N-1 

K=MODC Ks 2) +1 

READC2'K) PTS 

READCS? K)-- LCE 

CALL HAMFILCBSTR) 

CALL WINDOWCICEs 1, 2000) 
CALL DISPLACICEs 1s 2000s 25 1) 
CALL DISPLACPTSs 1s 2000s 25 2) 


CK KOK OK OK OR RK IO OK IK II OI OK OK IO KO A I IK I I KK IK 2k 2K 2k 2k ak ak ok ak 2k 


C 


PS THEMINDERLINE OK. ? 


CK OK RO OO ROK IK AK KOK IOI II IO I IOI I I I I I I I I a I kK a kk ok akc ak ok ak 2k 2k ak akc ok 


6 


107 


6 1 


CALL BGNPICC 3) 

WRITEC8s 107) 
FORMAPGIHIs” OK? Y OR N °) 
CALL ENDPIC 

READ(&s, 101) V 

CALL: DELPLE GL) 

CALL DELPIGGe) 

CALL DELPICC 3) 
CALIEWDELPLEGI ) 
IFCV-EQ-‘'N‘) GOTO & 

DO 61 K=1s, 2000 
PCEGKO =! CE CK) -PTSCK) 

CALL WINDOWCICEs 1s 200) 
CALL DISPLACICEs Is 2000, 151) 


C25 2K 2K KK OK 2 EK KK 2K OK OK 2K FE 2 KK 2K 2K 2 2c IK > 2K I 2K 2 2c 2 2 2K og 2 2K 2c 2c ig 2 2K 2g ok ok oc 2 KC 2K ig og 2K 2c is 2k 2 ic oi 2k 2c ok ok 2K ok ok 


C 


ES THE SRESUL BeGOOD ENOUGH» TO’ ‘GO? ON" TAPE ? 


CK 2K 2K 2K AK 2 KOK OK 26K 2K eo 2 KOK 2K 2K KK 2K oe 2K 2K 2c ig 2 iC 2 2K 2 2 2K 2g 2K 2c aK 2 KC ok 2K 2K 2 2g 2K 2 I 2K 2K 2K KK 2K 2K 2 2c 2K KEK ok 2K 2k 


103 


10 


CALL BGNPICC 3) 

WRITEC&s 107) 

CALL ENDPIC 

READ(8.101) V 

CALL DELPLCC 13 

CALE DEL PIG C3) 

CALL DELPIC(C99) 

IFCYeEQ-'N') GOTO 8 

CALL MTWRTC1,1ICEs 2000, IERR) 
M=M+1 

WRITEC1"° 1) ANCIL 

IFCIEOF-NE-U) GOTO 11 

CALL BGNPICC 3) 

WRITEC8>s 103) 
FORMATCIH1s*APREST? CY OR N)*) 
CALL ENDPIC 

READ(85101) V 
TFOVseNE«'Y*-AND-V-NE-~'N*) GOTO 10 
CALBSRELPIT CC 32 

IFCV-EQ-'N') GOTO 3 
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CORO OR OR OR OR IOI IO RK OI III OI AK IO IK I I I II I ICI I AK 2k ak 2K 2K 2k 2k a 2k 
C CLOSE OUTPUT TAPE AND REWIND BOTH BEFORE STOPPING 
COO ORO OI OK I IK I I IOI OK Ia I IK I a IK I a I a a a aK a aK I 2k 2k ke 2k ak 28 2k ak ak aK akc ok 2 
Ht IFCIEOF-NEeU) WRITEC 65 104) 
104 FORMATC' END OF DATA‘) 

CALL MTWTM(C 1) 

CALE Mien MiGels) 

CALL MTRWNDC 1) 

CALL MTRWNDC 0) 

STOP 
CORR ROR IO IRI KK IK III I OK I I ak aK IK kk ak 3k ak aK aK aI ak ak ake ak ak ak ak ak 2k 
G NOT GOOD ENOUGH. GO BACK AND DO IT AGAIN 
CORO OO I II RIO IK I I KR I IOI I OK OK I I a CI Kk ak kk 2k 2k ok ok akc ok ake 2k ak 
8 N=N-1 

GOTO 21 


END 
C2 2 2K 2K 2 2 I OK 26K 2k ic 2 2K 2K 2g 2K is ic 2 Ig 2c ic 2K ik 2K ik 2K 2g ak ig 2g 24 6 2c ig 2c 26 aK IC 2k 2k 2g 2K 2g 2 2c Ig 2 IE 2K 2c 2K 2K a 2K 2K 2K 2K 2g 2K 2K OK 2K ok 3K 


(CK 2K 2K 2K 26 2K 2 KK 26 2K > 2g 2K 2 2K 2 fk aK 2c 2k 2K Kg IK ok ic ik ik 2c 2k oe ic ic 2c ic ig 2c ak aK ik 2c ic 2c 2k 4c 3K 3k 2g 2 2c 2K 2K 2k 2k ofc 2 eK a oe oe 2K 2K 2k ok ok 2K 
C SUBROUTINE FOR READING A BLOCK FROM THE INPUT TAPE 
CK 2K 2K i aK 2 2K 2K 2K 2K IC 2 oC 2K oi aC 2 2K > 2K 2K 2K 2k a 2K 2 2 IC KK 2K 2k 2K 2K feo Ig ik 2K 2k 2c og ig 2c oi aK 2 2c 2k 2g 2K 2K 2c 2k 2K 2c 2g 2K 2K 2K 2k oe ok 2k 
SUBROUTINE TAPINCICEs Ns I EOF) 
INTEGER ICEC 1) 
DATA NB/1I/ 
IFCN-EQ-NB) GOTO 1 
K=N-NEB 
CALL MTSKIPCUs Ks IERRs I EOF) 
IFCIERR-EQ-0) GOTO 11 
WRITEC6, 100) N 
100 FORMATC' ERROR ON SKIPPING TO BLOCK's15) 
STOP 
11 IFCIEOFeNEeU) RETURN 
NB=N+ 1 
CALL MTREADC Us ICEs 2000, IERRs I EOF) 
DO 2 T=], 2000 
2 IPCECIJ=-I CEC?) 
IFCIERR-E@-e-0) RETURN 
WRITEC6, 101) N 
101 FORMATC* ERROR ON READING BLOCK',I15) 
STOP 
END 
CK KK OK KK KK OK RE 2K KKK KKK 2 2K 2 2K 2 KK OK 2k 2K 2k ic 2k 2k ok ak kc ke oo oo ok KK kK 2k oie 2k ak oC ocak 2K of Ko akc ake ae 
C38 2K IK KK 2K KK KK KK > 2K 2K 2K ok ok ok 2k 2k ok KK OK ok akc ok 3 2 2 KK 2K KKK ak a a 2k ak KK aK KK ok 3K ak ok ok ok ok ok 
C SUBROUTINE TO CORRECT DATA FOR PHASE SHIFTS AND NOISE SPIKES 
C2 2K KK 2K 2K 2K 2 2K 2 2K KK OK 2 OK 2K 2K 2k KK OK ko 2K oC ak ok KK 2 2K 2K 3k ak 2 ak ok 2k aK 2k kK a ak kK aK aK 2 ak oko ak ok ok ok ok ok ok ok ok 
SUBROUTINE CORREKCICEs OFSETsN) 
COMMON/A/W( 4) /B/TCC 2) 
DIMENSION XC 2)5YC 2) 
INTEGER LEDCe) aT CEC 1). OF SETs LsH 
DATA IED/2s,-1/7 
KA=OFSET 
CALL WINDOWCICEs 1s, 2000) 
CALL DISPLACICEs 1s 200U, 13 1) 
45 DO 44 1=131999 
J=ICECI+¢1)-I CECI) 
IFC J ¢LTe GUN ANDeJ-GTe-200) GOTO 44 
TCC 1D =CI-WC1))/¢CWC 2) -WC 1))d 
TCC 2) =WC 3) 
CALL WRI TOLC 2015 0s Us TCs 2) 
K=INWAITC-1l1-esI EDs IDAs; 1) 
LRGRDA SNE eG Of Oma 
JJ=I+1 
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DO 23 K=dJs eu0uu 


23 PCRCK) =1CECK)=J 
OF SET=OFSET<-dJ 
44 CONTINUE 


CALL BGNPIC(C 3) 
WRITEC8, 101) N 

Buri FORMATCIHIs’' TRY G R X. OR O's 10X51 4) 
CALE TENDRLG 

1 I=INWAITC-1-s, IEDs IDAs!) 
PECTIA.NE. *G!’) GOTO 2 
CALL REATOL(C 2015 05 0. TCs 2) 
XC 1) =TCC1)*CWC2)-WO1))4+W0 1) 
YC1)2=TCC2) * CWC 4) -W0 3) 04+WC 3) 
KI=IFIXCXC€1)) 
I=INWAIT(C-1-e5,I1ED) 
CALL REATOL(C 201505 0s TCs 2) 
XC Q)=TCC1I*CWC 2) -WCO1))4+WC0 1) 
YC2Q)=TCC 2) * CWO 4) -WC 3))+WC 3) 
Ke=CFIXCX€2)) 
PVeTRIXRCY Cl) =<¥ C2) ) 
DO 11 KK=K1s,Ke2 

1] ICECKK) =ICECK1) 
K2=Ke+1 
DO 12 KK=K2, 2000 

le ICECKK)J=ICEC KK +1V 
OF SET=OFSET+I1IV 
GOTO! 

2 TECIDA-NE.«*R* SAND. IDA-NE-*X'-AND-IDA-NE-'0") GOTO +1 
Pen DAsNEe 2.) GOTO 22 
CALL DEL RL CC 3) 
CALL TAPINCICE,N,s I1EOF) 


OF SET=KA 

DO 2! KK=15 2000 
el ICECKK) =I CEC KK) +0OF SET 
ea CALL DELPICCYY ) 


CALL WINDOWCICEs 1, 2000) 

CALL DELPICC 1) 

CALL DISPLACICEs 1, 2000s 14 1) 

IFCIDA-sEQ-'X') GOTO 45 

IFC IDA-EQ-'0O') RETURN 

GOTO ! 

END 
CK RR ROR OR ROR RO ROK OIRO IK GI IORI OR OK a koi kok ak 
CAR ORK OR ROK OO IO IORI IR IK IOI I I KK IK IK a ak aK ok 
C SUBROUTINE TO UNDERLINE THE DATA 
CRRA OO OR I RR OI IGRI IO IK I OI I OK a II I ak aC ok a ak aK ak 

SUBROUTINE UNDERLCPTSsLPT) 

COMMON /A/WC( 4) /B/TCC 2) 

PVVEGER OLE DCS) ePTSC1)s 1 TMPS bo) 

DATA ILED/2s-1/MAX/48/ 

CALL DELPIC¢ 3) 

CALL BGNPICC 3) 

WRITEC8, 100) 
10a FORMAT(1H1.s'UNDERLINE*) 

CALL ENDPIC 

CALL BGNPIC(C2) 

CALL ENDPIC 


6 ITMPC I) =IFIXCWC1)) 
LIMECS)=LPT 
K=3 
CALL DRAWCITMP.sK) 

] I=INWAIT(C-1-es1TEDs,1DAs 1) 
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LECT Dave PGi LR. Ge Oued 
LFiCLDAatlEe Gi) GOTO: 72 
LFCKeLTsMAX). GOTO iL 1 
CALL sDELPICG3) 
CALL .EBGNPICC 3) 
WRLTECS se t0)) 
map| PORMAT GIN Va. DELETE SOM Eve) 
CALE EN DELC 
GOTO. 13 
bad CALL REATOL CEU )s Ws WoT Cs.2) 
LTIMPCK) =LFIXC TCC 1) * CWC 2) -WO 1D) 4+WC 1) 
ITMPCK+ LO =1FIXCTCCE)* CWC 4) -WCG)) tWOCSD 
LE CLIMP CK). GT -1IMPCK=.2 il arG Ge Omae 
CALLE BDELPLECC 3) 
CALL) BGNPLGCS) 
WRITEGSs 102) 
Ie FORMATCIH1]s,'MOYVE ALONG‘) 
CALL ENDPIC 


GOrOn] 

Wee K=K+e 
CALLS DRAW GiaiM Pak 

es Gow Gaal 

2 PEG UC ACN Evers.) iG OO ae 
PeGKheE We oG.O Orel 
K=K-2 
CALL DRAWCITMP.sK) 

21 COP Os! 

a IEG VM IDYAG Ns VEY) Coe) SS 


GALL (REATOL Geu Tots. Oe Care) 
LPT=IFIXCTCC2)* (WC A) -WC 3) 4+WC 3D) 
GOLO26 
£5) CALL DELEL CCS) 
GALL BGNPLGC3) 
WRITECSoe1 U3) 
ANS FORMATCIH!s* TRY Bolsa Ds OR WO) 
CALL ENDPIC 
GOTrOe | 
4 CAL US VEEP LCGIS 
CALL. DELPICC2) 
CALE DEL PICe 3) 
CALE VEPPIECcg9) 
TCC 1)=0- 
TCC QD=U-5 
CALLAWRITOLC COU Ls As Us TCs 2) 
ITMPCK) =2000 
DIMPCK+ t= LIMPCK= 1) 
K=K-e 
DO 41 IT=tsKs2 
IX=I TMPCI+2) -ITMPCI) 
IY=I TMPC1+3)-1 TMPCI+ 1) 
Li = IMP Gl 
Le=LIMPGE+2) 
DOV a lea Li alee 
ITA=L-ITMPCI) 
At PTSCLIRELITMPCI+¢1) +1 F ITXC FLOATCIY) *FLCATCIAY/FUCATCInDD 
LPT=PTS¢ 2000) 
RETURN 
END 
LOR OOK ROR ORO ROR OI ORI IOI III aa a ake ke ak 
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CORIO COICO ICO OR IORI III I IK aK ak ak aK ak ak 
C SUBROUTINE TO DRAW PART OF AN UNDERLINE 
COICO IO IOI ISI a II IK aC kk a kK ak: 
SUBROUTINE DRAWCITMP.sK) 
COMMON /A/WC4)/B/TCC 2) 
INTEGER ITMP(C1) 
CALL DELPICC2) 
CALL BGNPICC 2) 
CALL\ LINICITMPC 1)5ITMPC2).0) 
CALL LINIRCUs Us 4) 
N=K-2 
IFCKeEQe-3) GOTO 2 
DO 1 I=35N32 
] CALL LINICITMPCI)sITMPCI+1)) 
CALL ENDPIC 
TCC 1) =C€C FLOATCITMPCN) )-WC1))/0WC2)-WC1)) 
TCC 2) =C FLOATCITMPCN+1))-WC3)0)/0WC 4) -WC 3) ) 
CALL WRITOLC 2015 05 Us TCs 2) 
RETURN 
END 
CK 2K 2K KOK OK 2K OK 2K CK OK 2 2g iG 2 2 2 kc 2 CK 2 2K KK 2K 2K 2K 2k ig oko 2 2g a 2K ic ke ok oa 2k ak aK 2k aI a ak ok 2c 2c ok 2 2c ok 2K 2k ok kc 
CK 2K 2K 2K 2 OK 2 2K 2K KK 2k 6 2K 2K 2K 2K 2K 2K 2K 2 2K 2K 2k 2k 2k 2 2k 2K ig 2K 2K oo 2 aK 2K 2g 2 2c i 2K 2k 2c oi 2c 2K a 2c oi kK 2 ok 2c ok 2 ok oie ok 
C SUBROUTINE TO DISPLAY DATA 
C6 KK 2K 2K 2K OK OK OK OK KE EK OK 2K 2K OK EK 2K 2K KK 2 i 2K oc ik 2c 2g 2 2S 2 2g 2K 2k 2g KK 2K kK oc ok 2k 2 ic oi 2k KK 2g KK ok 2K ok 2K 2k ok ok 2K 
SUBROUTINE DISPLACIAsLsHsKBYsNPIC) 
DIMENSION IAC!) 
INTEGER LsH 
CALL. BGNPICCNPIC) 
CALL LINICLsIACL) 5 0) 
DO 1 I=LsHsKBY 
1 CALL LINT<Is IACI) ) 
CALL ENDPIC 
RETURN 
END 
C5 2K 2 RK 2 2 2 2 2 2 2 C2 Eg 2c I Ck 2 i I 2 i 2 2k 2 a I I kK IS Ic 2 iC aig ke 2 ic akc 2K akc 2 fe ak 2 2c 
(CK KOK ARR OK OOK OK OK OOK IO OK OK OK OK IO OK OK IK IK IK OK IK a I I IK 3 2 kK I a IK ak ak 2k 2k ok akc 2k 2k ok 2k 2k 2k 2k 
C SUBROUTINE TO WINDOW THE DATA 
Ce 2 OK 2 2 A ER 2 2 A I A I 2 I 2 I 2 OK kk 
SUBROUTINE WINDOWCIA,L;H) 
COMMON /A/WC4) 
INTEGER IAC1)sYLsYHsLsH 
DATA IVAL/200/ 
YL=IACL) 
YH=IACL) 
DONT K=Loh 
TRePACK I «cLliterL2 YL=IACK) 
IFCIACK) eGT-YH) YH=IACK) 
1 CONTINUE 
WCO1)=L 
WC 2) =H 
WC 3) =YL-20U0- 
WC 4)=YH 
CALL WINDWCW) 
CALL BGNPIC(99) 
CALL LINECWC1).WC3)50) 
CALL LINIR(CUs,IVAL) 
CALL LINECWC1)505 0) 
CALL LINECWC2),0) 
CALL ENDPIC 
RETURN 
END 
2K KK 2 2 2K 2 2 2K 2K 2K eK 2 2K 2 2K 2k 2 kK 2 2k 2 2 2k i 2 2c oie ag oi 2 2K 2c 2g I aK 2K akc Ko 2 2k oe 2k 2s ic oi ok oc ako 2K 2k 2c ok 
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CO ORO I IK IO I I IOI IR I II II II IOI I CK GK I I IC IK IC aK ak a ak ac ak 2k akc a 
G SUBROUTINE TO TABULATE COEFFICIENTS FOR THE FILTER 
CORK OO IR OK OK III OI IK OK GK IK IO I I I OK I a 2k 2k IK IC I 3 2K 2k a 2K 2k 2 ak a 2k ke 2 akc 2 2c 
SUBROUTINE HAMTAB 
COMMON /C/TABL(C 200) 
DATA P/3-1415926536/ 
A=-20U.6 
DO 1 K=1s3,200 
TABL (CK) =Uc5« 54406 46* COSC P*A/S2QUUe ) 
l A=A+2 
RETURN 
END 
COO OO IORI IORI II I OK IR IORI IG OK I IO IIa Ia I a a I 2K 2k 2k ok 2 ake 2k akc ak ak 2k 
CoO OR OKI GO OK IGRI IOI I OK IK IO I I I I kK I 2K 2 2k 2 2k 2k 2k a 2 2k ke 2 akc ak ai 
C START RAMP FOR FILTER 
CORR KOR IO OIG OK ROOK OR IO III RII a a I I ak a 3 2K 25 2 5 25 2k 2k ak 2k ak 2k ak 2k ak 2k akc ak 
SUBROUTINE HAMBGN(CBSTR) 
INTEGER BSTR(C1) 
DO 1 I=13,100 
AI=I 
1 ESTRCI)=BSTRCI01) 
RETURN 
END 
CROOK IO ORO IORI RII IG RIG II OK a a I ak aK I ak aK 2k ak ak ak ak ak ak ak 
CO OOK OOOO IG OR I I IIR OK OI I a IK I kk a ak ak 2k ak ak a5- 
C END RAMP FOR FILTER 
CoRR KK OOK OK KK OK OOK IO I OK II I I IK I a OK OK ak ak 2k ok ak ak ok akc ak ak akc ak ak af 
SUBROUTINE HAMEND(C ESTR) 
- INTEGER ESTRC1) 
DO 1 IT=1]5,100 
ESTRCI+1)9=€C 100--FLOATCID *ESTRCOII/S1I00.~ 
l CONTINUE 
RETURN 
END 
CR IO IO OI ORR OR OC I IG II II I aK OK kK 2k aK 2k 2 2 2k a 2k aK akc a ak ak 
Co RIO I OR IGOR IOI III IK KR a IK IK I aK I I a IK Kak ak a ak ak ak ak ak ak ak ok 
€ COPY ROUTINE” FOR FILTER 
CoRR OK OK OK OR IK IO IO I I I I I I 2K I OK 2K I IK kk OK ak ak a aK ok akc ak ak ok ak ok 
SUBROUTINE HAMCOPCFPTSs ESTR) > 
INTEGER PTSC1)3ESTRC 1) 
DOe bh beds 100 
1 ESTRCI+1)=PTSCI) 
RETURN 
END 


CARR OR RRO OR ICO CR OR II OR I OR OR OI I IK CI ak a ak a ak 2 ak ak ak ak 
Ck AOC RO CI CIO IC IO OI IORI OR II I I OK IK kK GK kK kK a ak a ak ok ak ak 
C SUEFOUTINE ‘TO FILTER THE UNDERLINEs, IN PRACTICE@THIS (WAS 
C REPLACED BY AN EQUIVALENT ASSEMBLEF ROUTINE. 
CRI OOK I IC OI IO IC IO OK OK IORI ICR IOI RR OK I KK a aK a ak 2k ak 
SUBROUTINE HAMFILGBSTR) 
COMMON /C/TABL( 200) 
INTEGER BSTRC 1) 
DO i IT=l,e000 
SUM=0 
DO 2 K=15,200 
2 | SUM=SUM+BSTRCI+K=- 1) *TABLCK) 
1 BSTRCI)=SUM/108.- 
DO 3 T=1,2100 
3 BSTRC2ZE0)<I)=BSTRCeloi=1) 
RETURN 
END 
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DISPLAY PROFILE 


THERE 
A LARGE 
DISCONTINUITY 
BETWEEN 


REMOVE 
DISCONTINUITY 


NO 


REMOVE DISCONTINUITY 
DEFINED BY LIGHT PEN 


REDRAW PROFILE WITH 
ALTERATIONS SHOWN 


ALTERATIONS INCORRECT 
SO REDRAW ORIGINAL 
PROFILE 


REDRAW PROFILE 
WITH ALTERATIONS 
SHOWN 


RETURN 


Figure 10. CORREK subroutine 
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WAIT FOR 
COMMAND 


i ; athe 
"DELETE SOME" 
SEGMENTS? 


INSERT LINE 
SEGMENT FROM 
PREVIOUS POINT 
TO LIGHT PEN 


LIGHT PEN TO RIGHT 
OF PREVIOUS POINT? 


WRITE 
"MOVE ALONG" 


DELETE PREVIOUS 
LINE SEGMENT 


PLACE INITIAL POINT 
AT LIGHT PEN 
YES 


NO 


ERASE DISPLAY 


| 


| RETURN] 


Figure 11. UNDERL subroutine. 
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